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The magnetization dynamics in diluted magnetic semiconductor heterostructures based on 
(Zn,Mn)Se and (Cd,Mn)Te has been studied experimentally by optical methods and simulated 
numerically. In the samples with nonhomogeneous magnetic ion distribution this dynamics is con- 
tributed by spin-lattice relaxation and spin diffusion in the Mn spin system. The spin diffusion 
coefficient of 7 ■ 10~* cm^/s has been evaluated for Zno.99Mno.01Se from comparison of experimental 
and numerical results. Calculations of the giant Zeeman splitting of the exciton states and the 
magnetization dynamics in the ordered alloys and parabolic quantum wells fabricated by the digital 
growth technique show perfect agreement with the experimental data. In both structure types the 
spin diffusion has an essential contribution to the magnetization dynamics. 
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I. INTRODUCTION 

Heterostructures based on II- VI diluted magnetic semi- 
conductors (DMS), like Cdi_a;Mn^Te or Zni^^jMn^j^Se, 
are treated nowadays as model structures for spintronic 
devicesii^. A large variety of optical spectroscopy tech- 
niques has been conveniently used for investigation of 
their properties due to the strong optical transitions in 
the vicinity of the band gap being visible in, e.g., ab- 
sorption, reflection and emission. DMS materials are 
known for giant magneto-optical effects originating from 
the strong exchange interaction of free carriers with the 
localized magnetic moments of Mn^+ ionSi2,ii. 

The spin dynamics of the free carriers in DMS is con- 
trolled by the spin and energy transfer in the coupled 
system of the carriers, magnetic ions and phonon bath 
(lattice)^"—. The fast exchange scattering of the carriers 
on the magnetic ions (10~^^ — 10~^^ s) provides an effi- 
cient transfer into the Mn system, but the further transfer 
from the Mn spin system to the phonon bath is deceler- 
ated by the relatively slow spin-lattice relaxation (SLR) 
rate in materials with low Mn concentrations of a few 
percent. The control of the dynamical properties of the 
Mn spins, described by the SLR time, is essential because 
they may turn out to be a bottleneck for high-speed ap- 
plications. Therefore, fundamental knowledge about the 
dynamical properties of DMS and particularly about the 
dynamics of the Mn spin system is of importance. 

In n-VI DMS the SLR time of the localized Mn spins, 
TsLRj is controlled by concentration dependent exchange 
interactions between the Mn^+ ions. It varies by several 
orders of magnitude from milliseconds down to nanosec- 
onds with increasing Mn concentration x from 0.004 to 
0.11, see Fig. [T{c)^i^. Such a strong dependence is ex- 
plained by the role of Mn clusters as relaxation centers^ 



in the absence of SLR for isolated Mn ions. With increas- 
ing Mn concentration the number of clusters grows and 
the cluster structure evolutes in size involving more and 
more spins. This shortens a characteristic length for spin 
transfer from a single Mn ion to a cluster and also accel- 
erates the spin relaxation processes inside the clusters. 
It has been shown that the digital growth techniques of 
ordered DMS alloys allow for the same Mn concentra- 
tion a redistribution of cluster statistics, thus increasing 
the role of the Mn spin clusteraiS. The cluster formation 
results in a considerable enhancement of the SLR con- 
firming the leading role of clusters in the magnetization 
relaxation dynamics. 

There are some experimental evidences that spin diffu- 
sion in the Mn spin system plays an important role in the 
magnetization dynamics of II- VI DMS. Firstly, the strong 
dependence of the SLR time on the Mn concentration can 
be explained by a nonlocal model of the SLR in DMS 
where the spin diffusion from single Mn^+ ions to Mn 
clusters provides fast SLRii. Secondly, direct experimen- 
tal proof of the spin diffusion contribution to the magne- 
tization relaxation in heteromagnetic DMS structures has 
been reported for Zno.99Mno.01Se/Beo.93Mno.07Te struc- 
tures with different Mn concentration in the adjacent lay- 
ersi^. The slow SLR channel in the Zno.99Mno.01Se layer 
has been bypassed by the spin diffusion through the Mn 
spin system to the Beo.93Mno.07Te layer with much faster 
relaxation rate. Neither experimental nor theoretical re- 
ports allowing a quantitative evaluation of the spin dif- 
fusion parameters in DMS have been published, yet. 

The goal of this paper is to provide a detailed ex- 
perimental study of the spin diffusion in the Mn spin 
system of II- VI DMS and to evaluate the spin dif- 
fusion coefficient. We took advantages offered by 
Zni_a:Mna:Se/Bei_j,Mnj,Te heteromagnetic structures 
for a clear experimental emphasis of the spin diffusion 



contribution to magnetization relaxation dynamics. We 
developed a numerical approach to describe the magne- 
tization dynamics in structures with a modulated pro- 
file of Mn concentration and approved its validity for 
CdTe/Cdi_j;Mna;Te digital superlattices and parabolic 
quantum wells. 



II. EXPERIMENTALS 

Time-resolved magneto-optical studies of the mag- 
netization dynamics have been performed for various 
Zni_a;Mna;Se and Cdi-^^Mn^^Te based heterostructures 
grown by molecular-beam epitaxy on (lOO)-oriented 
GaAs substrates. Details of the structure design are 
given in the corresponding Sections in combination with 
the presentation of experimental data. In this paper we 
include results for: (i) three Zni_2;Mn2;Se/Bei_yMnj,Te 
multiple quantum well (MQW) structures with a type- 
II band alignment, (ii) three CdTe/Cdi^ajMn^-Te digi- 
tal superlattices (SL), and (iii) two CdTe/Cdi_j;Mna;Te 
parabohc quantum wells (QW). The CdTe/Cdi_j,Mn^Te 
structures have been fabricated by the digital growth 
technique^. 

We used an all-optical experimental technique devel- 
oped to measure the magnetization dynamics-^'^. It ex- 
ploits an internal thermometer of the Mn spin tempera- 
ture Tmii, which is provided by the high sensitivity of 
the giant Zeeman splitting of excitons (band states) to 
the polarization of the magnetic Mn ions in an external 
magnetic field. Taking into account that the changes in 
the spectral position of the photoluminescence (PL) line 
A£'pl(S) are proportional to the magnetization changes 
AM{B), we obtain by means of a magneto-optical tech- 
nique direct information about magnetization dynamics 
in DMS heterostructures. For that purpose the studied 
DMS structures have been exposed to an external mag- 
netic field, a pulsed laser has been used for heating the 
Mn spin system and the dynamical shift of the PL line 
has been measured by a time-resolved technique. 

The relaxation process of the heated Mn spin system 
to the bath temperature is described by the magneti- 
zation relaxation time Tmr which is contributed by the 
spin-lattice relaxation with the characteristic time Tslr 
and the spin diffusion. The relaxation times are equal 
(tmr = Tslr) for disordered DMS alloys and/or in the 
case when the spatial inhomogeneity in the Mn heat- 
ing does not contribute to the magnetization relaxation. 
These conditions are suggested to apply implicitly or ex- 
plicitly to most of the experiments, but this is not the 
case of our study, since the structures have been designed 
to evolve the spin diffusion contribution. 

Measurements were performed at low temperature 
T = 1.8 K for samples immersed in pumped liquid he- 
lium and in external magnetic fields S up to 10 T. The 
magnetic field generated by a superconducting split-coil 
solenoid was applied parallel to the structure growth axis 
(z-axis). The photoluminescence signal was detected 



for right-hand (cr^) circular polarization, which is the 
strongest component in applied external magnetic fields. 
A pump-probe technique with time-resolved detection 
of the PL signal excited by the probe laser pulses has 
been used. A pulsed YAG laser with a photon energy of 
3.49 eV (wavelength of 355 nm) was used as pump to pro- 
vide heating of the Mn spin system. It generated pulses 
with 5 ns duration at a repetition frequency ranging from 
100 Hz up to 5 kHz. Excitation power densities up to 
> 200 kW/cm^ were reached. A semiconductor laser 
(photon energy 3.06 eV, wavelength 405 nm) with exter- 
nal amplitude modulation operated as the probe laser. 
Laser pulses with a duration from 0.5 up to 100 |as were 
synchronized with the pump pulses. The time delay be- 
tween the pump and probe pulses was varied by a digital 
delay generator in the range from up to 10 ms. The 
excitation density of the probe laser has been kept be- 
low 0.1 W/cm^ in order to minimize its contribution to 
the heating of the Mn spin system. The time-resolved 
PL spectra were recorded by means of a gated charge- 
coupled-device (CCD) camera with a temporal resolution 
of < 2 ns. Avoidance of a spatially inhomogeneous exci- 
tation has been realized by overlapped pump and probe 
excitation spots with diameters of ~ 500 |J.m and by col- 
lection of the PL signal from their central part with a 
diameter of ^ 50 |J.m. 



III. MAGNETIZATION DYNAMICS IN TYPE-II 
DMS HETEROSTRUCTURES 

The DMS heterostructures Zni_2,Mn2,Se/Bei_yMnj,Te 
have a type-II band alignment with a potential en- 
ergy minimum for conduction band electrons in the 
Zni_:rMna:Se layers and for valence band holes in the 
Bei_yMnj,Te layers. The large difference in the band 
gaps, which are about 2.8 eV for Zni_2;Mna;Se and 4.5 eV 
for Bei_yMnyTe, and huge band offsets of > 2 and 
^ 1 eV for conduction and valence bands, respectively, 
prevent the penetration of the electron and hole wave 
functions into the neighboring layers. Details of the 
growth, band structure, and optical properties of these 
structures can be found in previous publicationsiiii^. 

We have studied three MQW samples with 
a type-II band alignment: sample ^1 with 

Zno.99Mno.01Se/BeTe (20/10 nm), sample #2 with 
Zno.99Mno.01Se/Beo.93Mno.07Te (20/10 nm), and sample 
#3 with Zno.99Mno.01Se/Beo.93Mno.07Te (10/5 nm). 
Each sample has ten periods with a respective layer 
thickness of 20, 10 or 5 nm. Sample #1 does not contain 
Mn in the BeTe layers. It serves as a reference structure 
to the sample ^2, which only differs by the presence of 
Mn in the Beo.93Mno.07Te layers. The main difference 
between the samples #2 and #3 is the thickness of the 
Zno.99Mno.01Se layers of 20 and 10 nm, respectively. 
This results in different distances for spin diffusion 
from the center of these layers to the interfaces with 
Beo.93Mno.07Te. 
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(a") BeTe ;c = 0.01 BcTe 
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FIG. 1: (a) and (b) Schematic band diagrams of 
Zni_j;Mnj;Se/Bei_yMni,Te quantum wells with a type-II 
band alignment. Dashed arrows show spatially direct (DT) 
and indirect (IT) optical transitions. Presence of the Mn ions 
in the Bei-yMn^Te layers, see panel (b), opens a new chan- 
nel for magnetization relaxation in the Zno.99Mno.01Se layer 
via spin diffusion into the Bei-yMn^Te layers, see text, (c) 
Spin-lattice relaxation time as a function of the Mn content 
in II- VI DMS structures based on Cdi-2,Mna;Te (closed cir- 
cles) and Zni-ajMn^Se (open circles) (Ref. @); the solid line 
represents the interpolation of this dependence. 



In Fig.[T][a,b) the band diagrams for the studied type-II 
structures are shown schematically. The spatially direct 
(DT) and indirect (IT) optical transitions are marked by 
arrows. In the photoluminescence spectra they appear 
as lines at ~ 2.8 and ~ 1.8 eV, respectively, as can be 
seen in the insets of Fig. [2j In our experimental con- 
ditions the photons, generated by the pump as well as 
the probe laser, are absorbed by the Zni_a;Mn2,Se lay- 
ers only. Hence, the indirect emission line is attributed 
to the holes scattered from the Zni_a;Mna;Se into the 
Bei-yMuyTe layers after the photogeneration. This hole 
relaxation process is fast and depends strongly on the 
width of the Zni_2:Mn2:Se layer— li^. As a result, the re- 
lative intensities of the direct and indirect emission lines 
in the samples #2 and #3 (pay attention to the multi- 
plication factor for the DT line of the sample #3) differ 
considerably. Since we are interested in the magnetiza- 
tion dynamics of the Zni_a;Mna;Se layers, we will focus 
on the giant Zeeman shift of the DT line in the following. 

Due to the strong exchange interaction of the free car- 
riers with the localized spins of the Mn^"*" ions the ex- 
citon transitions in DMS show a giant Zeeman shift in 
external magnetic fields. This shift equals to one half 
of the giant Zeeman splitting of the exciton commonly 
described by the modified Brillouin function, for details 
see, e.g., Eqs. (1) and (2) in Ref. [q]. The experimental 
results for the giant Zeeman shift in the three studied 
samples are depicted in Fig. [5] All three samples show 
very similar shifts saturating at the level of ~18-20 meV, 
which confirms that the static magnetic properties of the 
Zno.g9Mno.01Se layers are not affected by the presence 
or absence of the Mn ions in adjacent layers. This is a 
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FIG. 2: Giant Zeeman shift of the spatially direct PL line 
in different Zni_a;Mna;Se/Bei_t,Mnj,Te QWs. Experimental 
data are shown by symbols, the solid line is a fit for the sample 
#1 with parameters given in the text. In the insets the pho- 
toluminescence spectra of the samples #2 and #3 are shown. 



predictable behavior, as the Mn-carrier exchange inter- 
action is short-range and the carrier wave functions are 
strongly localized in the Zno.99Mno.01Se layersiiiS. The 
solid line in Fig. [5] is a fit by the modified Brillouin func- 
tion of the experimental data for the sample ^1. It has 
been done with the following parameters: ScS — 2.26, 
To = 1 K and x — 0.011. Similar fits performed for the 
samples #2 and #3 show that the Mn concentrations in 
the Zni_2;Mna;Se layers of all three samples differ by less 
than 8%. Therefore, their dynamical characteristics can 
be directly compared with each other. 

The dynamics of the magnetization relaxation after the 
pump pulse heating are illustrated for all three samples 
in Figs. [3] and |H In the sample #1 the relaxation time 
Tmr ~ 200 [IS coincides well with the SLR time known 
for Zno.99Mno.01Se QWs^. Indeed, for this sample no 
contribution of spin diffusion to the magnetization relax- 
ation is expected, since only the Zno.99Mno.01Se layers 
contain Mn ions. 

In comparison with the sample ^^1 a very strong re- 
duction of Tmr down to 6 |JS is found in the sample 
#2. Obviously, this is due to the spin diffusion from 
the Zno.99Mno.01Se layers to the Beo.93Mno.07Te layers, 
where the spin-lattice relaxation time is expected to be 
very fast (< 200 ns), see Fig. [TJc). Note, that very 
close results have been observed in similar samples in ex- 
periments, where pulses of nonequilibrium phonons have 
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FIG. 3: Temporal evolution of the changes in the spectral 
position of the PL line AEpl in Zno.99Mno.01Se/BeTe (20/10 
nm) MQW, sample #1. 



been used for the heating of the Mn spin systenJ^. 

A further shortening of Tmr down to 1.7 \is in the 
sample #3 by a factor of about 4 (in comparison with the 
sample #2) originates from a twice smaller distance for 
spin diffusion between the centers of the Zno.99Mno.01Se 
layers and the interfaces to Bco.93Mno.07Te. We conclude 
from these experimental data that the spin diffusion in 
Zno.99Mno.01Se can provide a spin transfer over distances 
exceeding 10 nm. We will show in the following Section 
IV that the results of Figs. [3] and |4] are sufficient for 
the evaluation of the spin diffusion coefficient. It makes 
heteromagnetic structures very suitable for quantitative 
studies of the spin diffusion in DMS. 



IV. NUMERICAL SIMULATIONS OF 
MAGNETIZATION DYNAMICS 

The magnetic relaxation time Tmr, extracted from the 
dynamical shift of the PL line, can be treated as the spin- 
lattice relaxation time Tslr only in case of uniformly 
distributed disordered magnetic ions. As it has been dis- 
cussed in Sec. Ill, the spin diffusion in the structures 
with a nonuniform Mn distribution may strongly accele- 
rate the magnetic relaxation. Generally, when the Mn 
concentration x{z) is a function of the growth direction 
z, the magnetization M{z) and the Mn spin temperature 
Tuniz) are also dependent on z during the magnetiza- 
tion relaxation process. Hence, the spectral position of 
the PL line, i?pL, is determined by the magnetization in 
the regions averaged by the wavefunctions of the carriers 
involved in the recombination. In this case the infor- 
mation about the temporal changes of the spatial profile 
of Tivin and, thus, about the local magnetization is very 
important for the adequate description of the experimen- 
tally obtained dynamics of A_EpL(t). We have developed 
a numerical model suitable for a direct comparison with 
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Zno.99Mno.01Se/Beo.93Mno.07Te structures with different 
layer thicknesses: 20/10 nm (#2 - open circles) and 10/5 nm 
(#3 - closed circles). 



experimental results in different samples. 



A. Concentration dependence of SLR dynamics 

The experimental results for the SLR dynamics of the 
Mn spin system are summarized in Fig. [Ijc) , where the 
SLR time is plotted as a function of the Mn concentra- 
tion. The open circles are data for Zni_a;Mna;Se-based 
structures, the data for Cdi_a:Mn3:Te are illustrated by 
the closed symbols^. The SLR times cover five orders 
of magnitude from 10~^ down to 10~® seconds as the 
Mn concentration varies from 0.004 up to 0.11. Such 
a strong dependence evidences that Mn-Mn interactions 
play a key role in the SLR dynamics. The results for 
both materials closely coincide demonstrating that the 
SLR of the Mn'^^ ions is rather insensitive to the ion 
host material. 

We turn now to the physical mechanisms responsible 
for the strong dependence of the SLR time 'rsLR(a;) on the 
Mn concentration. It is known that an isolated Mn^"*" ion 
does not couple with the lattice in a perfect II- VI semi- 
conductor crystal!^ and, therefore, the SLR time should 
be extremely long. This expectation can be explained by 
the following reasons: The electric field of the phonons 
does not act on the magnetic moments of the 3d electrons 
and also the magnetic field, induced by the temporal vari- 
ation of the electric field, is relativistically small. In prin- 
ciple, gradients of the electric field may interact with the 
quadrupolar moment of the magnetic ion which, in turn, 
interacts with the spin. But for the Mn^"'" ions the mag- 
netic d-shell is half-filled and its quadrupolar moment is 
equal to zero. In this case the dominating mechanism 
of the spin-phonon interaction is caused by the phonon 
modulation of the spin-spin interaction between neigh- 



boring Mn ions. The phonons modulate the distance be- 
tween these ions, hence, the local exchange field of an ion 
fluctuates. It is known as the Waller mechanisriJ^, which 
is obviously strongly dependent on the concentration of 
magnetic ions. Only the anisotropic exchange interaction 
is relevant because its operator does not commute with 
the Zeeman one. 

The spin-spin interactions for the Mn^+ ions are pro- 
vided by three mechanisms: the exchange interaction, the 
Dzyaloshinsky-Moriya interaction^SiSi and the magneto- 
dipole interaction. All of them decrease with increasing 
distance between the Mn ions, but the magneto-dipole 
interaction has the longest range. Therefore, despite be- 
ing relatively weak, it plays the dominant role for SLR 
dynamics in the limit of small Mn concentrations. The 
exchange and Dzyaloshinsky-Moriya interactions gain in 
importance at higher Mn contents, when the formation 
of Mn clusters is more probable. Both the number of 
Mn spins coupled in clusters as well as the typical cluster 
size progressively increase with growing Mn concentra- 
tion^^. This is the reason for the strong dependence of 
the SLR time on the Mn concentratioit^ i^^'^'^ . Within the 
framework of a simple model which accounts for the spin 
diffusion from Mn ions to Mn clusters, where the spin 
has an efficient relaxation, it can be shown that the SLR 
rate 1/tslr has a strong dependence on the Mn content 
with a power law between x^ and x^ |ll]. This is in good 
qualitative agreement with the experimental results of 
Fig.IHc). 

According to Fig. [Ijc), one can conclude that for fixed 
external conditions the Mn concentration x mainly deter- 
mines the SLR dynamics. This establishes a relationship 
between the static and dynamic magnetization, e.g., be- 
tween the giant Zeeman splitting value and the SLR time 
- both being dependent on the Mn content. This restrict- 
ing relationship can be relaxed or even canceled by the 
growth of heterostructures with modulated Mn content 
and/or with modulation doping by donors or acceptors 
providing free carriers. In Section V few examples of such 
structures are regarded, where the magnetization dynam- 
ics has been considerably accelerated without changing 
the static magnetization. 



B. Model 



the antiferromagnetic Mn-Mn interactions contributing 
to the overall magnetizatio n^^ i ^^ . In the equilibrium con- 
dition the Mn spin temperature is equal to the tempera- 
ture of the phonon bath: Tyin = ^l. 

The model should take into account the processes 
of spin diffusion and spin-lattice relaxation in the Mn 
spin system. It should finally describe the relaxation of 
the magnetization changes |AAf (ATmh)! — \M{TMn) — 
M(Tl)| induced by an external pertubation in a constant 
external magnetic field. In this non-equilibrium situa- 
tion the spin temperature difference is given by ATmh = 
Tun — Ti^- In DMS with a homogeneous Mn concentra- 
tion X, which does not depend on spatial coordinates, the 
magnetization dynamics in the one-dimensional case can 
be described by the following equation: 



d_ 

di 



|AM|=ifdiff(x)^|AAf| 






|AM| 
TsLR(a;) 



(2) 



Here Kiig{x) is the spin diffusion coefficient. In the case 
of a dipole-dipole interaction between the Mn^+ ion spins 
the dependence of K^is on the Mn concentration x can 
be expressed as: 



-firdiff(a;) ^ Kqx 



,4/3 



(3) 



where Kq is the spin diffusion constant, see Eq. (3.25) in 
Ref. [13. 

In consideration of Eq. ([1]), the Eq. ^ can be written 
in terms of ATmii^ 



d (9^ 



|AAf| r dMjTMn) 

tSLR(2:) V ^^Mn 



._ (4) 

In the case of small deviations from the equilibrium state 
and/or if Af is a linear function of Tmhi the last term be- 
comes equal to ATmh/tslr- It should be noted that for 
the ranges of B and Tmh, regarded in the present pub- 
lication, the relation between the magnetization changes 
AAf and ATun is close to linear. Therefore, we will sug- 
gest that the relaxation of ATmh is simply determined 
by TsLR. 

The equation ^ can be generalized for a nonuniform 
Mn distribution along the growth direction z of the hete- 
rostructurc x{z): 



The magnetization M{B,TMn) in DMS is a function 
of the magnetic field strength and the Mn spin tempera- 
ture^: 



M{B,TMn) = 5^Ba;A'^0'S'eff(x)B5/2 



5gnBB 



2kB{TMn + Toix)) 



(1) 

Here Nq is the inverse unit-cell volume, B5/2 is the modi- 
fied Brillouin function, g — 2.01 is f;-factor of the Mn^+ 
ions, /cb is the Boltzmann constant and /ib is the Bohr 
magneton. Scs{x) is the effective spin and Tq(x) is the 
effective temperature. These two parameters account for 



d_ 
dt' 



AT, 



Mn 



^- (i^diff(z)^-ArMn 

oz \ oz 



AT, 



Mn 



tslr(^;) 



(5) 



In the model we suggest that the spin diffusion coefficient 
Kdisiz) is determined by the local Mn concentration x{z) 
according to Eq. ([3]). Additional assumption is that the 
SLR time is determined by the local Mn concentration 
x{z). For the dependence TsLR(a;) we used the interpo- 
lation function shown in Fig. [ijc) by the solid line. 

In order to compare the experimental and numerical 
results the calculations have been performed in two steps. 
Firstly, we have solved numerically the one-dimensional 



partial differential Eq. ([5]), which describes the dynamics 
of the spatial distribution of T^n in the sample with a 
nonuniform concentration of magnetic ions, as an initial- 
boundary value problem. We have chosen the uniform in 
space distribution of Tmu as the initial condition. The so- 
lution gives us information about Tmh in different points 
in space and at different time moments. An example of 
such kind of calculation, which will be discussed in detail 
below, is shown in Fig. |6l 

The second step comprises the numerical solution of 
the Schrodinger equations for electrons and holes. The 
potential profiles for carriers in DMS depend on the ex- 
ternal magnetic field strength and the Mn spin tempera- 
ture due to the giant Zeeman effect. The corresponding 
energy shifts of conduction and valence bands are de- 
scribed by: 



AE,^h)iB,TMn) 



(X, 



(h) 



mg 



MiB,TMn) 



(6) 



where iVgae and NqOCh are the exchange constants for 
the conduction and valence band. In Zni-^jMn^jSe they 
are equal to 0.26 and -1.31 eV, and in Cdi-ajMn^^Te they 
are equal to 0.22 and -0.88 eV, respectively^. The Mn 
spin temperature distribution TMn{z,t), obtained as nu- 
merical solution of the Eq. ([5]), allows us to determine 
the exact potential profiles for electrons and holes in an 
external magnetic field at different times after the pulsed 
heating of the Mn spin system. For every moment of 
time we solved numerically the Schrodinger equations to 
obtain the electron and hole energies. This, in turn, al- 
lows us to calculate the changes in the spectral position 
of the PL line and to compare it with the experimental 
dataof A£'pL(i)- 



Numerical results for type-II DMS 
heterostructures 



In this Section we present the results of numeri- 
cal simulations of the magnetic relaxation dynamics in 
the type-II Zno.99Mno.01Se/Beo.93Mno.07Te (20/10 nm) 
MQW (sample #2). The experimental results for this 
structure are shown in Fig.|4l The Mn profile x{z) in this 
structure, which we used to solve Eq. ([S]), is a step-like 
function with the values of 0.01 and 0.07 in the adjacent 
layers. The SLR times Tslr in Zno.99Mno.01Se (~ 200 ^s) 
and Beo.93Mno.07Te (< 200 ns) differ by more than three 
orders of magnitude. In that case the Beo.93Mno.07Te 
layers serve as a very efficient SLR channel for spin ex- 
citations which diffuse from the Zno.99Mno.01Se layers 
to them. The results of the numerical simulations be- 
come nearly independent on the spin diffusion inside the 
Beo.93Mno.07Te layers and are predominantly controlled 
by K(iis in the Zno.99Mno.01Se layers. Therefore, this 
type of heteromagnetic heterostructure is very reliable 
for an evaluation of ii'diff. 

In Figure \5\ the calculated dependence of Tmr on Kdiff 
is shown by the thick solid line. The magnetic relaxation 
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FIG. 5: Calculated dependence of the magnetic relax- 
ation time Tmr on the diffusion coefficient Kdift in the 
Zno.g9Mno.01Se layer of the Zno.99Mno.01Se/Beo.93Mno.07Te 
(20/10 nm) QW. In the limiting cases of small and large 
i^diff the magnetic relaxation time is determined by tslr in 
Zno.99Mno.01Se and Beo.93Mno.07Te, respectively. The dashed 
line shows the experimental value of Tmr ~ 6 (xs. 



times for this dependence were used as characteristic 
times for the calculated A£^pl (t) dynamics. These times 
were taken at the point where the normalized A£'pL(t) 
is reduced to 1/e = 0.37 of its original value. Thus, 
the results in Fig. [3] can be directly compared with the 
experiment. As expected, the extremal values of K^is 



a^^^3 







FIG. 6: Calculated temporal and spatial distribution of the 
Mn spin temperature in the Zno.99Mno.01Se/Beo.93Mno.07Te 
QW (20/10 nm). Initially, Tivin = 10 K is uniform in 
space, the final temperature is Tmu = Tl = 2 K, K^is = 
7 ■ 10~* cm^/s in the Zno.99Mno.01Se layer. Thick solid lines 
show the temperature behavior close to the interfaces. 
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FIG. 7: (Color online) Calculated spatial distribution of the 
Mn spin temperature in the Zno.99Mno.01Se/Beo.93Mno.07Te 
QW (20/10 nm) at different time delays after the pump pulse. 
Panels (a) and (b) correspond to different values of Kdm in 
the Zno.99Mno.01Se layer. The initial and final temperatures 
are 10 and 2 K, respectively. 



— ^ and -^ 00 reflect the identity Tmr = Tslr, in the 
Zno.99Mno.01Se and Beo.93Mno.07Te layers, respectively. 
The dashed arrow in Fig. [5] indicates the experimental 
value for Tmr « 6 |as, measured for the sample #2 (see 
Fig. |4|). From its crossing point with the solid line we 
evaluate the spin diffusion coefficient in Zno.g9Mno.01Se 
with rather high accuracy to i^diff ~ (7±1.5)-10~* cm^/s. 
Then, according to Eq. ([3]) the diffusion constant is 
Kq k, 3.5 • 10~^ cm^/s. To the best of our knowledge, 
this is the first experimental evaluation of the spin diffu- 
sion parameters in DMS. 

Let us make now a closer look at the calculated spatial 
and temporal profiles of Tmh in the heteromagnetic struc- 
ture #2 (see Fig.[S]). The initial condition for this numeri- 
cal problem is that the Mn spin system of the whole sam- 
ple is heated homogeneously up to Tmh = 10 K, which 
can be seen as a fiat temperature level at short time de- 
lays. In our experiments such heating is provided by the 
interaction of the Mn ions with photogcnerated carriers. 
The subsequent magnetic relaxation in the Mn spin sys- 
tem is provided by the SLR and by the spin diffusion. 
The latter is responsible for the spatially inhomogeneous 
cooling of the Mn spin system in the Zno.99Mno.01Se 
layer, which is clearly visible in Fig. [Sj At time delays 
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FIG. 8: Temporal evolution of the Mn spin temperature in 
the center of the Zno.99Mno.01Se layer calculated for different 
values of the diffusion coefficient A'diff. The initial and final 
temperatures are 10 and 2 K, respectively. 



exceeding 20 |as all Mn spins reach the equilibrium tem- 
perature of the phonon bath, which has been chosen to 
Tl = 2 K. 

For more detailed insight the spatial profiles of Tyi^{z) 
at different time delays are shown in Fig. El^a). In fact, 
these profiles are cross sectional areas of Fig. [6] at fixed 
time delays ranging from 0.1 up to 100 \xs. In order 
to demonstrate the importance of the spin diffusion for 
the magnetic relaxation in this structure we present in 
Fig. [T^b) the calculated results for the case, where the 
spin diffusion is absent, i.e., for i^diff = 0. It is consistent 
with an independent relaxation of Tmh in the adjacent 
layers with the corresponding SLR times. The tempera- 
ture profiles in each layer are flat. Also, the relaxation 
of Tmii for, e.g., Mn ions located in the vicinity of the 
center of the Zno.99Mno.01Se layer is considerably slower 
than in the case of a flnite diffusion coefficient, as shown 
inFig.ElJa). 

Figure[5]dcmonstratcs the dynamics of Tmh in the cen- 
ter of the Zno.99Mno.01Se layer which is strongly depen- 
dent on fCdiff in this layer. For fCdiff = the spin diffusion 
is absent and the relaxation process is solely controlled 
by the SLR time in the Zno.99Mno.01Se layer, which ap- 
proximately amounts to 200 |j,s. On the contrary, if the 
value of Xdiff — >■ 00, and, thus, the spin diffusion time 
does not restrict the magnetic relaxation time, the Mn 
spin temperature drops very fast to the bath value of 
Tmii = Tl = 2 K, being only controlled by the SLR time 
(< 200 ns) in the Beo.93Mno.07Te layers. The calculated 
dynamics of Tmh with iiTdiff = 7 ■ 10~* cm^/s, which cor- 
responds to the experimental situation of the sample ^2, 
is shown by the dashed line in Fig. [51 



V. MAGNETIZATION DYNAMICS IN 

CdTe/CdMnTe STRUCTURES ENGINEERED BY 

DIGITAL GROWTH TECHNIQUE 

The digital growth technique has achieved about 
atomic precision by using molecular-beam epitaxy. It 
has been implemented for GaAs/AlAs heterostructures, 
where during the growth the wide-band-gap AlAs lay- 
ers with a certain thickness (frequently a submonolayer) 
have been introduced into narrow-band-gap GaAs layers 
at strictly predefined positions^S,. This technique allows 
a highly variable engineering of the shape of the confining 
potential in quantum well structures. Also, by growing 
periodic structures with layer thicknesses of the order of 
a single monolayer a digital alloy (DA) can be fabricated. 
Its properties can significantly differ from the properties 
of the common disordered alloys, e.g., Al2;Gai_a;As. 

The digital growth technique has been successfully 
established for DMS heterostructuresi^i^. For the 
CdTe/Cdi_j;Mn2;Te material system digital alloys and 
parabolic quantum wells have been grow n^°i'^° . In these 
structures the band gap profile is controlled by the imple- 
mentation of the Mn ions, and, therefore, the magnetic 
properties become also modulated. It has been shown 
experimentally that the partial smoothing of the digi- 
tally grown Mn profile considerably reduces the possible 
modifications of the static magnetic properties due to 
the Mn intermixing with neighboring CdTe layers^iiS^. 
Typically, the latter are measured optically via the giant 
Zeeman splitting of the excitons. However, the magneti- 
zation dynamics can be changed drastically^*^. Such a dif- 
ference in the static and dynamical magnetization prop- 
erties can be explained by different mechanisms. The 
static magnetization is mainly contributed by the para- 
magnetic Mn spins, while the magnetization dynamics 
critically depends on the concentration of the Mn clus- 
ters which serve as spin relaxation centers. 

In this Section we analyze the Mn spin dynamics in 
two types of CdTe/Cdi_a;Mna;Te heterostructures: digi- 
tal alloys and parabolic QWs. We will demonstrate that 
the spin diffusion in the Mn spin system plays a key role 
in its dynamics, which is also simulated numerically. 



A. CdTe/CdMnTe digital alloys 

Let us analyze the experimental results for three 
samples with CdTe/Cdi_a;Mn2:Te digital alloys re- 
ported in Ref. [13] ■ In these samples 23-nm-thick 
Cdo.95Mno.05Te/CdTe DA layers (corresponding to 
about 70 monolayers) are confined between nonmagnetic 
Cdo.8Mgo.2Te barriers. According to the thickness of the 
magnetic and nonmagnetic layers, given in the monolay- 
ers, the structures are labeled in the following way: the 
first sample 1 x 3DA consists of 18 periods of 1ML/3ML, 
the second sample 2 x 6DA has 9 periods of 2ML/6ML, 
and the third sample 3 x 9DA has 6 periods of 3ML/9ML. 
The choice of the layer thickness provides the same av- 
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FIG. 9: Giant Zeeman shift (a) and dynamical shift of the 
exciton PL lines AEpi^{t) during the magnetic relaxation (b) 
in three digital alloy samples. The sample design is shown 
schematically in the inset. The experimental data in both 
panels, shown by symbols, are taken from Ref. [10]. Solid lines 
are results of numerical calculations. For the modeling the 
values of the diffusion coefficient KdiB ~ 3.5 • 10~^a;^' ''[cm /s] 
and Mn intermixing length lo — 0.3 nm were used. The 
dashed line in panel (b) represents the calculated results for 
the 2 x 6DA sample, neglecting spin diffusion (Kdia ~ 0). 



erage Mn concentration xda ~ 0.013 in the three DA 
samples. 

The experimental data for the excitonic giant Zeeman 
shifts and magnetic relaxation dynamics, measured in 
Ref. [lOl, are shown by the symbols in the Figs.|9l[a) and 
[9jb), respectively. One can see that the Zeeman shifts 
of all three samples follow closely each other with rather 
weak differences in the saturation levels at high magnetic 
fields. This reflects the fact that the static magnetization 
is mainly controlled by the average Mn concentration in 
the DA. 

On the contrary, the magnetic relaxation dynamics of 
the digital alloys strongly depends on the details of the 
real (smoothed) Mn profiles. As one can see in Fig. [9Kb), 
the magnetic relaxation times shorten from 9.5 ps in 1 x 
3D A to 2.5 |xs in 3 X 9DA. They considerably deviate from 
the relaxation time of 27 \is, measured in the disordered 



Cdo.gssMno.oisTe reference alloj*iS. 

We have simulated numerically the magnetic relax- 
ation dynamics in these digital alloys. For that purpose, 
we smoothed the nominal technological profile by imple- 
menting the Mn intermixing, i.e., we took the convolution 
of this profile with an exponential broadening function'^^ : 



ff(^) 



exp(z/?o - 1), z <Iq 
0, z^lo 



(7) 



Then we solved the Eq. ([5]) by using the diffusion coeffi- 
cient in form of Eq. ^ with Kq — 3.5 • 10^*^ cm^/s, as it 
was estimated in Section IV. We realized that the Mn in- 
termixing length Iq strongly controls the relative ratios of 
the three Tmr values in different DA samples, which of- 
fers a high precision for the evaluation of Iq w 0.3 nin. It 
is important to stress here, that the magnetic relaxation 
in the digital alloys does not follow a single exponen- 
tial law, however, it can be rather closely fitted, compare 
solid lines and symbols in Fig. IHlJb). It is also interest- 
ing that, despite the very small layer thicknesses in the 
DA, the consideration of the spin diffusion in the Mn 
spin system was important. The dashed line in Fig. Hljb) 
shows the calculated result for the 2 x 6DA sample in 
the absence of spin diffusion {Kdiff = 0). It considerably 
deviates from the experimental data. 

Using the Mn intermixing length ^o = 0.3 nm, eval- 
uated from the best fit of the magnetic relaxation dy- 
namics, we have calculated the giant Zeeman shift and 
received good coincidence with the experimental data, as 
the solid lines and symbols in Fig. ^a) indicate. This 
makes us confident that the used numerical approach 
properly accounts for the main factors which control the 
static and dynamic magnetization in the studied digital 
alloys. 
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FIG. 10: (a) The thin line represents the technolog- 
ical profile of the Mn concentration in a 25.4-nm-thick 
CdTe/Cdi-2:Mna;Te parabolic QW. The effective averaged 
Mn concentration is given by the thick line (see text), (b) 
The efi'ective conduction and valence band profiles for the 
parabolic QW with Xmax = 0.68. Gray areas show the squared 
moduli of the wave functions of electrons and heavy-holes in 
the ground states. 



B. CdTe/CdMnTe parabolic QWs 

The digital growth technology enables the fabrication 
of structures with a continuously varying Mn concentra- 
tion and potential profile along the growth direction. One 
of the interesting examples of the technological capabil- 
ities is a parabolic CdTe/Cdi_a:Mna;Te QW, where the 
confinement potential for the carriers in the QW has a 
parabolic dependence on the coordinate along the growth 
direction. Details of the growth and magneto-optical 
properties for such QWs can be found in Refs. 30ll33ll34| . 

We have measured the magnetic relaxation dynamics 
for two parabolic QWs with different maximal Mn con- 
tents Xmax of 0.13 and 0.68. Both QWs have a thickness 
of 25.4 nm, which corresponds to 82 monolayers (one 
monolayer of CdTe is equal to 0.324 nm). The Mn shut- 
ter was continuously opened for the growth of the barri- 
ers and pulsed during the growth of the parabolic QW. 
The growth of each QW was divided into 41 steps. In 
each step the Mn cell was opened for the time needed to 
obtain a given value of Mn concentration, more details 
can be found in Ref. [SQl. The technological Mn profile 



averaged over each monolayer is shown by the thin solid 
line in Fig. llOf a) . After its strong smoothing one can re- 
ceive a parabolic profile represented by the thick line. In 
the real structures the Mn intermixing is not sufficiently 
strong to reach the monotonous Mn distribution shown 
by the solid line, it only smoothed the sharp boundaries 
in the nominal Mn distribution. 

The carrier wave functions, for both the electrons in 
the conduction band (CB) and the heavy-holes in the 
valence band (VB), penetrate through the digital lay- 
ers. As a result, their quantum well confinement poten- 
tials can be well approximated by the parabolic functions 
shown in Fig. llOf b) by thick lines. The square modulus 
of each carrier wave function, calculated in the real po- 
tential with Xmax = 0.68 and by considering the Mn in- 
termixing length ^0 = 0.3 nm, is shown by the gray area 
in Fig. [TUT b) . One can see that the electrons and holes 
at the lowest quantum confined energy levels are local- 
ized close to the center of the parabolic QW. The details 
on the energy spectrum and exciton giant Zeeman split- 
ting in CdTe/Cdi„a;Mn3;Te parabolic QWs are given in 
Refs. [Hill. 
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FIG. 11: Magnetic relajcation dynamics in 25.4-nni-thick 
CdTe/Cdi_i:Mna;Te parabolic quantum wells with different 
a^max- Solid lines are calculation results with A'diff = 3.5 ■ 
10~^a;^'''[cm /s] and Iq = 0.3 nm. The dashed line shows the 
calculated results for the sample with Xmax = 0.68 without 
taking into account the spin diffusion (A"diff = 0). 



The experimental results on the magnetic relaxation 
dynamics, measured from the dynamical shift of the 
exciton PL line AE'pl, are shown by the symbols in 
Fig. [TTJ The characteristic times Tmr for the samples 
with Xmax = 0.13 and 0.68 differ considerably from each 
other, being ^^ 1 |xs and ^ 100 ns, respectively. It can be 
explained by the role of the Mn clustering, which controls 
the TsLR, and by the Mn spin diffusion. 

One can expect that in these relatively thick structures 
with a gradual variation of the magnetic properties from 
the QW center to its edges the spin diffusion in the Mn 
spin system should play an important role for the mag- 
netic relaxation. In fact, its significance for the sample 
with Xmax — 0.68 is pointed out by comparing the re- 
sults of the numerical calculations with (the solid line 
in Fig. [TT|) and without (the dashed line) regarding the 
spin diffusion. The very long relaxation tail of the dashed 
line dependence (note the logarithmic time scale) is due 
to the long SLR times in the region of low Mn concen- 
tration near by the QW center. This region is important 
for the dynamics of the exciton giant Zeeman shift, as 
the electron and hole wave functions are primarily con- 
centrated on the QW center, see Fig. [TOlb). 

Furthermore, the numerical results, presented in 
Fig. [TTJ do coincide with the experimental data for the 
parabolic QWs with very different Xmax = 0.68 and 0.13 
without changing the calculation parameters (Kdis — 
3.5 • 10-^a;'*/3[cmVs], ^o = 0.3 nm). This fact is very 
promising that the magnetic relaxation dynamics in such 
complicated structures can be well described within the 
frame of a relatively simple model and without fitting pa- 
rameters. It also makes us confident that the suggested 
functional dependence of the diffusion coefficient on the 



Mn content, given by Eq. ([3]), is well suitable for the 
modeling of the magnetization dynamics. 



VI. DISCUSSION AND CONCLUSIONS 

Note, that even in bulk disordered DMS samples with 
a homogeneous distribution of Mn ions the spin diffusion 
can play a significant role for the magnetic relaxation 
dynamics. This fact can be attributed to at least two 
reasons: (i) Spatially inhomogeneous heating of the Mn 
spin systems, e.g., by a tightly focused laser beam. Based 
on typical parameters for the studied materials, we esti- 
mate that the scale of such an inhomogeneity should be 
smaller than 1 |j.in. (ii) Intrinsic instabilities of a highly 
excited Mn spin system, which may lead to the forma- 
tion of domains with "hot" and "cold" spin tempera- 
tures, whose size and dynamics will be governed by the 
spin diffusiort^"— . These conclusions are equally true 
for heterostructures containing layers with only one Mn 
concentration, e.g., DMS quantum wells confined by non- 
magnetic barriers. 

It is known, that the SLR dynamics of the Mn ions 
is strongly accelerated by an increase of the lattice 
temperatures^ and also becomes faster in stronger ex- 
ternal magnetic fields^^i^. Clearly, these effects would 
modify the magnetic relaxation dynamics in heteromag- 
netic heterostructures, which remarkably depends on the 
absolute values of tslr. 

In addition to that, the spin diffusion can be a dom- 
inating factor for the magnetic relaxation dynamics in 
epitaxially grown DMS quantum dots. Such quantum 
dots are typically grown on the top of the (about one 
monolayer thick) wetting layer, whose effective Mn con- 
centration (due to the very small thickness, see discussion 
in the Section for digital alloys) is considerably smaller 
than in the quantum dots. Experimental data available 
for such dots show a strongly nonexponential magnetic 
relaxation, which has been tentatively explained by the 
spin diffusion between the Mn spin systems of the quan- 
tum dots and the wetting layeri2ii£. DMS quantum dots 
have a faster SLR compared to the wetting layer, but 
they are characterized by a very limited capacity of the 
spin reservoir. In this case, the spin diffusion would con- 
trol the heating of the Mn spin system in contrast to 
the cooling scenarios, we have considered in the present 
paper. 

In conclusion, we have performed detailed experimen- 
tal and theoretical investigations of the magnetic relax- 
ation dynamics in different types of DMS quantum wells. 
We have shown that the spin diffusion in the Mn spin 
system plays an important role in this dynamics, which 
is especially evident for the structures with a nonuni- 
form Mn distribution. The diffusion coefficient i^diff = 
7 • 10~® cm^/s has been evaluated for Zno.99Mno.01Se. 
The results for the parabolic QWs suggest that the func- 
tional dependence of the diffusion coefficient on the Mn 
content in terms oi Kdis ~ 3.5- lO^^x*/'^ [cm'^/s] received 
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from the theoretical approach of Ref. [27| is suitable for 
the modeling in a wide range of Mn concentrations. 
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